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Abstract 

Scattering of light by the photonic clusters made of small particles 
is studied with the help of the quasiperiodicity condition and the local 
perturbation method. The analytical expression for the field scattered 
by the cluster is presented and the conditions of resonance scattering 
are found. The conditions of the zero scattering by the cluster are also 
found. 

1 Introduction 

Light scattering is a classical problem of electrodynamics which with emer- 
gency of photonic crystals became actual topic with many exiting applica- 
tions [I], [2]. The scattering by photonic crystals is difficult to treat analyti- 
cally because of they complex structure and large number of particles or cells 
involved in calculations. In practice, real photonic crystals are clusters with 
finite dimensions and finite number of particles (or cells). In many theoretical 



*The paper was submitted to JOSA B. 
'''e-mail: prosentsov@yahoo.com 



1 



works the clusters were approximated by infinite crystals and electromagnetic 
fields inside the clusters assumed to be quasiperiodic (see for example works 
[3], |1] and references therein). Even in experimental works the same technic 
was used successfully to explain observed results |5]. Alternatively, the light 
scattering by photonic clusters was studied rigorously by using the finite dif- 
ference time domain method [B] and the local perturbation method [7], [S] 
(often called also discrete dipole approximation or coupled dipole approx- 
imation method [9]). These numerical methods require number crunching 
machines with substantial computational power and memory. The principal 
limitation of the numerical methods is that final answer is a "naked" num- 
ber without a functional dependence. Some kind of trend can be obtained 
after several computations, functional dependencies however, can not be pre- 
dicted. In the local perturbation method, for example, the solution of the 
large system of linear equation is required and final result includes sums over 
all particles [7]. Analytical solution can be obtained in principle, however 
it will be not transparent for understanding. Under certain conditions (see 
for example work [8]), the cluster can be treated as infinite crystal where 
the quasiperiodicity exists. Due to the quasiperiodicity, the field should be 
calculated in one point only. This drastically simplifies calculations in the 
local perturbation method and the final results become more transparent. 
Firstly, it was done for one-dimensional infinite crystals [10] and afterwards 
it was generalized to three-dimensional structures [TT] . 

In our work we studied the light scattering from the spherical cluster 
by using both the local perturbation method and the quasiperiodicity. In 
contrast to previous works we obtained much simpler analytical expressions 
for the fields and analyzed resonance scattering. The resonance frequencies 
of the scattered field were calculated and resonance width was discussed. It 
was found that the field scattered by the spherical cluster has "cloaking" 
frequencies at which the field is absent. 
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2 The formalism 



The formalism we use is presented in many works (see for example [7]-[13]) 
and we only briefly present it here for convenience and consistency. Consider 
the photonic cluster made of identical particles which characteristic size L is 
small compared to incident wavelength Aq. The electric field E propagating 
in the host medium filled with the particles is described by the following 
equation [H] 



(A - V ® V + E(r) + k r^' V /(r - r„)E(r„) = S(r), (1) 
where 



^ I 27r CO ^ ^ I 1, inside particle 

Ao c ' I 0, outside particle 

Here A is Laplacian and V nabla operators, ® defines tensor product, fco is a 
wave number in the host medium {u is the angular frequency and c is the light 
velocity in vacuum), Esc and Eq are the permittivity of the particles and the 
medium respectively, / is the function describing the shape of the scatterers, 
and S is the field source. Note, that the equation ([1]) is an approximate one 
and this is true only when the small scatterers {k^L ^ 1) are considered. This 
equation is already solvable, however it will invoke solution of by system 
of linear equations for each frequency u. In addition, the final result will 
include summation over all particles that reduces solution transparency 
even further. The solution of the Eq. ([T]) can be significantly simplified when 
the field E is quasiperiodic, i.e. when 

E(r + r„) = E(r)e^''o'"'\ (2) 

We note that all particles are supposed to be identical for the condition ([2]) 
to be valid. By using relation ([2]) the equation ([1]) can be written in the 
following form 
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E(r) = Eo(r) + k^^-^^^^ (T+ E(0) e^'""^<^'n(r), (3) 



where 
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E(0) = E(ro)|r„=o> r„, = J^mJ'd,, (4) 

i=l 

Here / is the 3x3 unitary tensor in polarization space and r„ = Yl^=i fn^^i 
is the radius vector of the n-th particle, where are integers (m° = 0) and 
dj are the basis vectors of the crystal. The field E(0) is the field at the 0-th 
particle and / is the Fourier transform of the function /. The incident field 
Eq is created by the source S in the host medium and it is not important 
for our consideration (the expression for Eq is presented by formula ( 1A1I1 
in Appendix). The formula ([3]) is rather general one and it describes the 
field in the medium with photonic cluster of arbitrary form made of small 
particles of arbitrary form. The formula shows that the total field E is a sum 
of the incident field Eq propagating in the medium and the field due to the 
scattering by the cluster formed by particles. The interference between 
the scatterers as well as their resonance properties are taken into account by 
only one term E(0) that should be found by solving the system of 9 linear 
equations obtained by substituting r = tq = into (jS]). The formula ([3]) 
can be simplified even further when the distance between the observer and 
an n-th scatterer (-R„) is large, i.e. when i?„ ^ L. In this case the integral 
can be calculated approximately. We note also that the integral can 
be calculated exactly at least for the spherical particles. When Rn ^ L 
integration in Eq. ([3]) gives 
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where 



i?„ = |r-r„| > L. (7) 

Here i?„ is the distance between the observation point and n-th scatterer and 
V is the scatterer's volume. 

It is interesting to note that the expression describes the field inside 
infinite photonic crystal when N ^ oo. In addition, the formula explicitly 
shows that even small losses in the host medium (small positive imaginary 
part of ko) limit the number of interacting particles in the infinite photonic 
crystal. 

The formulae and (jH]) can be used with two requirements: the quasiperi- 
odicity and small size of the particles. The second requirement can be over- 
come by further subdivision of the particles. When the scatterers are not 
sufficiently small {koL > 0.2, for example), they can be subdivided into 
smaller particles with size Lgub such that k^Lsub < 0.1. 

It should be emphasized here that the main assumption of the work is the 
existence of the quasiperiodicity (expression ([2])) for a finite size photonic sys- 
tem. Obviously, the assumption is not correct for the lossless systems where 
light freely propagates from one boundary to another. Weak interference 
between the scatterers (due to small optical contrast or large distance, for 
example) and small losses allow to treat finite system as infinite one because 
of the effective absence of the boundaries. In this case the expression ([2]) 
should be correct. 

2.1 Fields in scalar approximation 

In this subsection the results obtained in scalar approximation are presented. 
These results are valid when light depolarization is negligible. We note that 
the results are important for general understanding and comparison with 
vector case. For definiteness we assume that the cluster is made of spherical 
particles with radius L (for particles with other shapes the similar results 
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can be obtained). In accordance with the formula the field in scalar 
approximation is 

E{r) = Eo{r) + klL' ^''\ £(0) e"^^''"-^, (8) 

Seo ^ 

where 

Em = ^ (9) 

and 



(r — } f,ikor„ 



n=l 

a = ^2 (gsc - ^o) (i^,2fcoL/3), r^=K\. (11) 

The expressions OH]) and f lTU]) already show that resonance properties of the 
field E{0) are specified by the parameters of the 0-th particle via a and by 
the rest of the cluster via the sum in Eq. ( !T0|) . The contribution of the 
rest of the cluster is controlled by its structure and varies from significant 
to negligible (as shown in Eqs. (|T2l) and (fl6ll below). The sums in the 
expressions ([8]) and ( 1T0|) can be calculated numerically, however the final 
result will be not transparent for analysis. Alternatively, these sums can be 
evaluated analytically by using the integration (see formulae flA2p and (IA3p 
in Appendix) when k^di <^ 1. This condition is extremely convenient for our 
consideration, because for close packed clusters (when di = 2L) it coincides 
with the assumption required for the method we use, i.e. koL ^ 1. For the 
cluster it is naturally to use the condition r ^ rn (observer is far from the 
cluster). The final expression for the field in the medium with the spherical 
cluster is 

Eir) = Eoir) + A;gL^/3 ^(0)^^, (12) 
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where 




C = MA + c^/2). 



(13) 



D 



l-a-2-K 



{Ssc - gp) L"^ i kl<fl^ - sm2(e) + ihd 
3eo d? I +i sm(20/2 - 2^ 



(14) 



A;o(A + c//2). 



(15) 



The expression (fT2|) is the field in the medium filled with small spherical 
particles forming spherical cluster of radius A. The particles are positioned 
in cubic lattice with period d. The formulae (fT2ll - (jnil is the main result of 
this subsection and we will discuss them thoroughly. 

Consider first the scattered field. The formulae f|T2|) and f|T3|) show that 
the scattered field oscillates as function of the scattering wave number ksc 
and the radius of the cluster A. The expressions ( |T2i) and ( |T3i) show also 
infiuence of the scattering vector k^c. When k^c = (backward scattering) 
P ^ N and the field is proportional to the cluster volume NL^. The formula 
(fT2|) also shows that the scattering field is zero when /3 = 0, i.e. when 



Taking into account that k^^d^ <^ 1, the solution of the equation (fT6|) can 
be found in the following form 



where M is a positive integer. The existence of the solution means that the 
cluster is actually invisible (cloaking effect) at the "cloaking" frequency Ud 
which is 




(16) 




(17) 
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We note here that the expression ( 1T8|) exphcitly shows that the cloaking 
frequencies of the photonic cluster depend on the cluster radius A and the 
direction of the scattering wave vector (defined by fj,). 

Consider now the resonance properties of the scattered field. The impor- 
tance of the Eq. (fT^ is that Re{D) = defines the resonance frequencies of 
the field scattered by the cluster. For the dielectric scatterers the equation 
for the resonance frequencies can be expressed in the form 

2er, ,9 9 / ttL\ AttL^ r,, ^ . . 

klL' l-^]+--^sm'{0. (19) 



(Bsc- Bo) "V QdJ 3 (P 

Note that when the period d of the cluster increases, the expression (fT9i) re- 
produces equation for the resonance frequencies of the single particle. Con- 
sider two most interesting cases. 

I. When sin(^) = the solution of the Eq. (IT^ is 

The formula (12U]) clearly shows that the resonance frequency of the single 
particle is the principal contribution and the period of the cluster only tunes 
the resonance frequency of the scattered field (keep in mind that d > 2L). 
In addition, the resonance frequency is defined from the condition sin(i^) = 
and it is 

ttM c 

where M is a positive integer. 

L 
d 



II. When f sin(0 > k'^d'^ the solution of the Eq. is 



arcsin ( ± 



3d3 £0 
2nL^ (esc-eo) 



B^ (A + d/2) 



(22) 
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The formula fl22p explicitly shows that the resonance frequency of the field 
scattered by the cluster is periodic and rather complex function of the cluster 
parameters (period d, scatterer size L, optical contrast Esc — £o) and the 
frequency is inversely proportional to the cluster size A. We note here that 
the resonance will occur when the following condition is satisfied 



3d^ 



< "° ^ ^ 1- (23) 

2ttL^ [Esc " ^o) 

The condition (123!) shows that for close packed cluster (when d = 2L), for 
example, the condition e^c > S^o should be satisfied. 

Finally, we consider the resonance width 7 which near the resonance 
(when CO = cjr) is defined in the following way 



\lm{D{u;))\ 



dRc{D) 




did 





(24) 



The resonance is sharp when 7(co'r) ^ c^r- In accordance with the formula 
( l24l) and (fT4ll the resonance width can be presented in the form 



7(w) 



kpL I 2L 
3 



sin(2g) 
2 



1 



(25) 



The expression 0251) explicitly shows that the resonance width is sensitive to 
the cluster size A ~ dN^^^, period d, and characteristic size of the particles L. 
The formula fl25|) also confirms that the resonance width of the field scattered 
by the cluster with extremely large period {k^Ld ^ ^N^^^) coincides with 
the resonance width of the field scattered by the single particle. When the 
cluster is relatively dense {d ~ 2L) and ^ sin (2^) > 1, the resonance width 
takes the form 



|sin(2e)| 

i.e. the resonance of such cluster is not sharp. 



(26) 
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2.2 Fields in vector case 



The vector case will be considered in a similar manner as the scalar one. 
Without the loss of generality we can assume that the wave vector of the 
incident wave kg is directed along z direction. This assumption simplifies 
the calculations of the sums in Eq. After summation in formula (EJ the 
field in the far zone [k^r ^ 1) can be presented in the following form 



E(r) = Eo(r) + "-^[l-X® l) E(0), (27) 

1 = r/r 

where the parameter /3 is given by the formula (fT3|l and the components of 
the field -E'j(O) are 

E,(0) = ^^^, (j = x,y,z). (28) 
Here the denominator Da is described by the following formulae 



D 



x,y 



1-X 



1-X 



(^.c-.o)4vrL3 / 41 (-1-1 + ^ + 4,) 



3£o 

{Esc - gp) 47rL3 
3£:o 



Hi 



X 



X x" / 



(29) 



x=kod/2 



\x x^ x" / 



x=^ 



x=kod/2 



When the dimension of the cluster is relatively large (^ ^ 1), the formulae 
(129|) can be simplified and written in the following form 



^0) / T ,„2 r 2 „• 2 , 3 r 3 



^ = (_3) J - ^oL' - t^KL' ]A = ko{A + d/2\ 



where 
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_ j ze/2-e2^V4 + 7/12, j = x,y 
'~\ .(l + e-^)/e-2A J = z ■ ^ ^ 

As well as in scalar case the scattered field oscillates with the cluster's size and 
vanishes when /3 = 0. In addition, the influence of the observation position 
is important via the scattering wave number kgc and the depolarization term 
/-I® I. 

The equations for the resonance frequencies are defined from the condition 
Re(Dj) = 0. The expressions fl30|l and fl3T|) show that the imaginary part 
of D^ y grows with the size of the cluster A preventing a sharp resonance. 
At the same time, the imaginary part of Dz decreases with the size of the 
cluster making the sharp resonance possible. Below we will study only the 
sharp resonance. The equation for the resonance frequencies is defined from 
the condition 'Re{Dz) = and the equation has the following form 

esc + 2eo ^2^2_8vrL3_ ^^^^ 



{Bsc - Bo) ' 3d^ 
In this case the resonance frequency is 



- ^ / Bsc + 2eo 8ttL^ 

The solution (133!) shows that the resonance frequency of the relatively large 
cluster (when the field is quasiperiodic) is actually modified resonance fre- 
quency of the single scatterer. When the period d of the cluster grows, the 
modification tends to zero very fast. Note that for the close packed clusters 
{d ~ 2L) the modification is essential and can not be neglected. 

It should be emphasized that the expression ( l33i) is very important and 
elegant one: it describes the resonance for the finite structure and does not 
contain the size of the structure. This effect can be explained by the usage 
of the quasiperiodic condition suitable for infinite structures. 
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3 Conclusions 

The analytical expressions for the field scattered by the spherical photonic 
cluster made of small spherical particles have been presented. The scattered 
field and its resonance properties have been studied. The resonance frequen- 
cies of the cluster were calculated. It has been shown that the resonance 
frequency of the scattered field is a function of the cluster's period and the 
characteristic size of the scatterers. 

The "cloaking" frequencies at which the cluster does not scatter have 
been found and they have been calculated explicitly. 
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5 Appendix 

5.1 The electric field created by a source 

The electric field Eg created by the source S in vacuum can be found by 
using the following formula (see for example jl5] ) 




(Al) 



where G{r, r') is the vacuum Green's tensor which is 




(Al.l) 
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5.2 The series summation via integration 

Consider the sum 

7V-1 

a=Y,e"''''\ {xn = nd), (A2) 

n=l 

where n is an integer and d is the period for particles placed at a;„. When 
kod <^ 1 (long wavelength approximation), the function e^''°^" varies slowly 
and the sum flA2p can be replaced by integration in the following way 



a = V e"'''^-^ ~ ^ / e'^'^^dx. (A3) 



n=l " 



Here we used the smallest space bin 5xn = d. 
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